Classical Phase Space Revealed by Coherent Light 
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We study the far field characteristics of oval-resonator laser diodes made of an AlGaAs/GaAs 
quantum well. The resonator shapes are various oval geometries, thereby probing chaotic and 
mixed classical dynamics. The far field pattern shows a pronounced fine structure that strongly 
depends on the cavity shape. Comparing the experimental data with ray-model simulations for a 
Fresnel billiard yields convincing agreement for all geometries and reveals the importance of the 
underlying classical phase space for the lasing characteristics. 
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Introduction — In the past two decades, quantum 
chaos has proven to be a successful concept in under- 
standing, characterizing, and predicting the behavior of 
mesoscopic systems [11,12]. Originally used to study the 
quantum analogue of classically chaotic hard- wall billiard 
model systems, it also explains the behavior of realistic 
systems of various shapes and character: The statistics 
of Coulomb-blockade peaks in quantum dots [2| , the ge- 
ometry dependence of the weak localization peak 0], or 
the level statistics of microwave billiards [4(. The stud- 
ies on quantum chaos have mainly focused on the quest 
for universality from the view-point of statistical physics. 
Therefore, despite the evident importance of the sys- 
tem's underlying classical phase space for the behavior 
of the quantum or wave mechanical analogue (based on 
the analogy between Schrodinger and Helmholtz equation 
its specific structure cannot be reconstructed from 
the traces it leaves in typical observables like energy level 
or wave function statistics. 

In this Letter we shall see that detailed information 
about the classical phase space can, however, be ex- 
tracted from the far field radiation characteristics of oval- 
microcavity laser diodes. The observed far field pattern 
(FFP) depends very sensitively on the shape, that is, 
on the system's underlying classical phase space, and we 
convincingly support this idea by numerical simulations. 
It is interesting to note that information about the clas- 
sical phase space is revealed by the coherent light ema- 
nating from the lasing microcavity. 

The oval-billiard family is a well-known model system 
with interesting properties [B|, @, @] in which each oval 
shape is characterized by a shape parameter 8 = 28' /L, 
cf. Fig. [lib). The curvature is, by construction [see 
Fig. [TJb) and, e.g., Ref. for instructions], discontinu- 
ous. This was found to be crucial for understanding the 
phase-space evolution upon variation of 8 Q. Increasing 
8, the phase space changes from fully chaotic (8 = 0, sta- 
dium shape), via mixed (0 < 8 < 1, various oval shapes) 
to integrable (5 = 1, disk). Figure [2] shows the phase 
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FIG. 1: (Color online) An oval microlaser with 8 = 
2<5'/L=0.45. (a) Scanning Electron Micrograph, (b) geom- 
etry, and (c) phase-space coordinates. 



space of the chaotic and integrable limit and the mixed 
phase spaces for 8 G [0.45,0.85] that consist of chaotic re- 
gions as well as regular islands [§] . The simultaneous ex- 
perimental and theoretical investigation of a 5-dependent 
mixed phase space has not, to the best of our knowledge, 
been done before and is at the core of the present study. 

To understand the dynamics of the system, let us con- 
sider a light ray hitting the resonator boundary at polar 
angle <f> under an angle of incidence %, cf. FigQJc), such 
that total internal reflection takes place (|sinx| > 1/n). 
In the circular cavity (8 = 1), the light ray will remain 
confined inside the cavity forever by means of conserva- 
tion of angular momentum, and form whispering gallery 
mode (WGM) tori. To a certain extent, this scenario will 
remain valid even if 8 is slightly decreased. The invariant 
WGM tori are then perturbed to Kolmogorov-Arnold- 
Moser tori, and the phase space becomes mixed. Note 
that for 8 >0.76 the chaotic regions are disconnected (the 
chaotic sea is split into two disjoined regions) whereas 
for 8 <0.76 a singly connected chaotic layer penetrates 
the phase space, bypassing the embedded regular islands 
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FIG. 2: (Color online) Poincare surface of section (SOS) il- 
lustrating the classical phase space of hard-wall oval billiards. 
The critical lines | sinxl = 1/n are shown by light- blue lines. 
In the top raw, the limiting cases of a chaotic (8 = 0, sta- 
dium) and an integrable (8 = 1, disk) phase space are shown; 
for 5=0.45, 0.76, 0.8, and 0.85 the phase space is mixed. To 
further visualize the differences in the phase-space structure, 
the evolution of 50 rays with random initial conditions in the 
interval < 4> < 0.1 and 0.28 < sinx < 0.32 is followed for 
25 reflections (superposed lighter points). 

H, [3, cf. Fig. m The smaller S the more chaotic the 

cavity becomes, reaching ergodicity in the limit S = 0. 

The above mentioned properties of ray-dynamics in 
oval billiards imply in particular that when S exceeds 
values of ~ 0.76, light rays are trapped in the dis- 
joined chaotic regions that lie above the critical line 
| sin x| = 1/n, unable to leave the cavity. Accordingly, 
the emitted light intensity is expected to decrease around 
6 ~ 0.76 because the lasing modes will be strongly con- 
fined inside the cavity and leakage will be evanescent. 

In order to elucidate if this property of classical phase 
space could be detected in the emitted light, we actually 
fabricated these oval-shaped semiconductor laser diodes 
for various values of 5, one of which is shown in Fig. 1(a). 

System — The resonators were made of metal organic 
vapor phase epitaxy grown graded index seperate con- 
finement AlGaAs/GaAs single-quantum well wafer with 
n=3.3 the effective index of refraction. The typical sys- 
tem length is L ~ 50 fim. Applying the reactive-ion- 
beam etching technique allows one to realize specific oval 
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FIG. 3: Total emitted intensity of lasing microcavities 
vs. shape parameter r5. Note the intensity drop (marked by ar- 
row) around 8 ~ 0.76 caused by the appearance of a transver- 
sal barrier in phase. 

cavities characterized by a shape parameter 5 with ex- 
tremely smooth and vertical boundaries, see Fig. 1(a); 
the surface roughness is less than 1/10 of the lasing wave- 
length A ~ 850 nm. Lasing operation of these devices was 
achieved at room temperature by pumping with a pulsed 
injection current of 500 ns width at a 1 kHz repetition 
rate. 

Experiments — In Fig. [3] we plot the total (polar-angle 
integrated) far field output power, measured at an injec- 
tion current of four times the threshold current, as func- 
tion of the shape parameter 6. The pronounced drop in 
the emitted light intensity around r5 < 0.76 (marked by 
arrow) is a direct consequence of the above-mentioned ap- 
pearance of a transversal barrier in classical phase space 
and the change of the lasing modes from refractive modes 
into evanescent, WGM-type modes. This dip is also a 
clear sign that not individual modes, but the underlying 
classicalphase space as a whole (representing a wealth of 
modes [8j) determines the behaviour of the system. 

In the following, we will shift our focus from the total 
output power to the angular distribution of the emitted 
radiation in the far field, i.e., the FFPs. The experimen- 
tal FFPs of oval lasing diodes with various 8 are shown 
in Fig. 2] (black curves). FFPs were measured in a dis- 
tance of d = 70 mm by scanning a photo-detector with 
a window of w = 11 mm width around the laser diode 
covering an angular range of 200° in </>ff with a resolution 
A0 ff = 2arctan(w/2ri) ~ 9°. 

The red curves in Fig. [4] are the results of numerical, 
ray-model based simulations (explained in more detail 
below) that convincingly support the experimental re- 
sults. A uniform background was substracted in all cases 
in order to exclude effects of spontaneous emission that 
is assumed to be isotropic; the (absolute) minimum of 
the experimental curves is set to zero by definition. 

Model — The theoretical model used here is based on 
the ray picture. This is motivated by the large size pa- 
rameter ~ 2imL/\ ~ 1000 of the experimental system; 
the ray model can be expected to yield reliable results 
in this regime, d EH A large number (250000) of rays 
is started with unit amplitude and random initial condi- 
tions covering the whole phase space. The dynamics of 
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FIG. 4: (Color online) Experimentally observed FFPs (far 
field intensity Ig vs. angle </>ff, black curves; same injection 
currents as in Fig. [3] were used) for various shape parameters 
8 (see insets). The ray-model simulations are superimposed 
(red curves); all data was averaged over a 9° window. The 
simulation data is offset by the average residual background 
of the experimental data, normalization is with respect to 
the peak heights. Apart from the side peaks, the agreement 
experiment-theory is convincing and the drastic changes in 
the FFP upon tiny changes in the shapes are nicely reflected 
in the simulation. 



each ray is governed by Snell's and Fresnel's law [llj that 
complement the condition of specular reflection. The far 
field radiation characteristics is reconstructed from rays 
in the steady phase-space distribution [H, [12] , an invari- 
ant object uniquely related to the underlying phase-space 
structure and characterized by an exponential decay of 
the internal light intensity vs. trajectory length. The ini- 
tial transient regime where the internal intensity decays 
faster than exponentially due to the presence of, e.g., 
bouncing ball orbits, is discarded [§]. That the linear 
ray-model works so well for describing non-linear, lasing 
light is related to the large size-parameter and the fact 
that not a single, but a multitude of modes is lasing. Q 

Results — The observed FFPs exhibit a clear and 
prominent minimum-maximum structure that is remark- 
ably well reproduced by the ray-model simulations for 
all 5, cf. Fig. [U Apart from the edges of the broad max- 
ima where side peaks are sometimes missing, the simu- 
lated FFPs capture the decreasing width and fine struc- 
ture of the central part, in particular its extreme shape- 
sensitivity, in a semi-quantitative manner. 

For all 5 considered here, the preferred emission direc- 
tion coincides with the shorter cavity axis (4>g ~ 90°). 
Note that this direction is shifted by 90° in compari- 
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FIG. 5: (Color online) From classical phase space to the FFPs 
of microcavities (ray-model based simulation) . The transmit- 
ted Fresnel intensity (color scale) is shown in the near field 
(lower panels, sine of angle of incidence \ vs - polar angle an- 
gle <f> at transmission) and far field (upper panels, resulting 
far field angle 4>fi vs - polar angle <j> of ray origin). The region 
of violated total internal reflection is enclosed by light lines in 
the lower panels. The intensity distribution in the lower pan- 
els corresponds well to the hard-wall results (cf. the lighter 
points in Fig. [2}; their translation into the FFP in the upper 
panels and in Fig. 2] is evident. The inset for 5 = 0.45 illus- 
trates how the details of the observed structure are related to 
shape-specific orbits as the one shown (see text). 



son with the stadium shaped microcavity laser § due to 
dynamical eclipsing [13]: For 5 > 0, regular (bouncing 
ball) islands around the stable fixed points sin x = and 
<f> = 0°, 180° impose constraints on all other orbits. Light 
rays will still emerge mainly from the regions of highest 
curvature, but nearly tangentially, i.e., in direction of the 
shorter cavity axis, cf. also Figs. [5] and [5l 

In the remaining part of this Letter, we will focus on 
the fine structure of the FFP maxima and show how the 
underlying classical phase space expresses itself in this 
signature, cf. Figs. 2] and [5] That the ray-model can ef- 
fectively reproduce the fine structure and width of the 
central part of the FFP maxima becomes especially evi- 
dent when looking at 6 = 0.76 where the side peaks are 
clearly separated and the plateau-like structure of the 
central part is well described, similar arguments apply to 
the wide maximum with the slight dip for S—0A5. The 
power of our simple model is also nicely illustrated when 
comparing the FFPs for S — 0.8 and 0.85 which are strik- 
ingly different despite the very similar cavity geometry: 
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For 6 — 0.8, a distinct minimum in the central plateau is 
evident in both experiment and theory which is missing 
for 6 — 0.85, again in both data sets. 

FFP and classical phase space — We now provide nu- 
merical evidence that these distinct differences in the fine 
structure of the FFP maxima originate in the structure 
of the underlying phase space. We already discussed the 
Poincare SOS of oval (hard wall) billiards in Fig. [2] and 
turn now to the properties of light leaving the open sys- 
tem, cf. Fig. 03 The Fresnel intensity transmitted into far 
held direction <pg, originating from a certain near held an- 
gle (f> and an angle of incidence x prior to transmission, 
is represented in color scale histograms for the near and 
far held intensity (lower and upper panels, respectively). 
To this end the respective spaces, (smx,4>) and (<j)ff,(j>), 
were divided into cells in which the Fresnel intensity of 
transmitted light rays was collected using the ray model 
described above. 

Superimposing the classical phase space from Fig. [2] 
on the light distribution in the near held (lower panels) 
reveals that for all S light rays leave indeed exclusively 
from the high curvature regions. Moreover, the larger 5, 
the closer to the critical lines | sin x\ = l/ n are the angles 
of incidence prior to transmission. That this is a charac- 
teristic property of the underlying classical phase space 
becomes even clearer when comparing the signatures in 
Fig. [5] with the fingerprint of light rays started at the 
high curvature region with near-critical incidence in the 
Poincare SOS, see the lighter points in Fig. [2] 

For 5=0.45, the structure in the emmiting regions is 
richer than Fig.[2]would suggest. This is related to the in- 
tricate interplay between geometry, trajectory, and Fres- 
nel's law and illustrated in the inset of Fig. [5j To this end 
50 rays (marked by white circles and numbered 0) leaving 
the cavity in the high intensity region close to the crit- 
ical line were followed for the next reflections (marked 
by numbers, colors and white trajectory are guides for 
the eye); those rays emit only little intensity at the first 
bounces. Intensity is again emitted refractively at the 
next subcritical reflections occuring 6 and 7 bounces later 
(arrows in inset). The Poincare SOS signature of these 
bounces is marked by green squares (the other numbers 
complete the Poincare SOS schematically): the emer- 
gence of the line-like structure in the near field Fresnel 
intensity (originating at bounce 7) becomes evident (sym- 
metry considerations complete the line pattern). Note 
that trajectories started in the intensity gaps couple out 
refractively after few bounces but not once the steady 
phase space distribution is reached. 

The upper panels of Fig. [5] show how the near field 
patterns are translated into the far fields, the correspon- 
dence is evident. Since the FFPs shown in Fig. |4] are 
obtained by summing the Fresnel intensities in the upper 
panels horizontally (over all <f>) , we conclude that the ob- 
served far-held characteristics of (even) the lasing cavity 
is intimately connected with, and can nicely be explained 



by, the classical phase-space structure. Moreover, origin 
and evolution of the fine structure of the FFP maxima 
in Fig. |4] are easily understood in the representation of 
Fig. El Each FFP maxima in Fig. |4] consists of two con- 
tributions corresponding to light leaving from the two 
regions of highest curvature. For 0.45 < 5 < 0.8, the two 
contributions are separated in tfiff causing a multipeak- 
structure of the FFPs, whereas for the only slightly larger 
shape parameter S — 0.85 they (almost) merge, in excel- 
lent agreement with experiment. 

Conclusion — We fabricated oval-shaped resonant mi- 
crocavity lasers of various shapes, and lasing operation 
was successfully established. We showed that the emit- 
ted coherent light uniquely reflects the characteristics of 
the underlying classical phase space. We employed a ray- 
based model that allows us to relate the observed peaks in 
the FFPs convincingly and unambiguously to the phase- 
space structure. In particular we find that the formation 
of a transversal barrier in the underlying classical phase 
space around S < 0.76 can be detected by coherent light. 
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